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ABSTRACT 

Polarization induced by cosmological scalar perturbations leads to a typical 
anisotropy pattern, which can best be analyzed in Fourier domain. This 
allows one to unambiguously distinguish cosmological signal of polarization 
from other foregrounds and systematics, as well as from polarization induced 
by non-scalar perturbations. The precision with which polarization and 
cross-correlation power spectra can be determined is limited by cosmic variance, 
noise and foreground residuals. Choice of estimator can significantly improve 
our capability of extracting cosmological signal and in the noise dominated 
limit the optimal power spectrum estimator reduces the variance by a factor of 
two compared to the simplest estimator. If foreground residuals are important 
then a different estimator can be used, which eliminates systematic effects 
from foregrounds so that no further foreground subtraction is needed. A 
particular combination of Stokes Q and U parameters vanishes for scalar 
induced polarization, thereby allowing an unambiguous determination of tensor 
modes. Theoretical predictions of polarization in standard models show that 
one typically expects a signal at the level of 5-10/iK on small angular scales and 
around 1/iK on large scales (/ < 200). Satellite missions should be able to reach 
sensitivities needed for an unambiguous detection of polarization, which would 
help to break the degeneracies in the determination of some of the cosmological 
parameters. 



Subject headings: polarization; methods: data analysis; cosmology: cosmic 
microwave background 
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1. Introduction 

Anisotropies in cosmic microwave background (CMB) are now widely accepted as the 
best probe of early universe, which can potentially provide information over a whole range 
of cosmological parameters (Jungman et al. 1996; Zaldarriaga, Spergel & Seljak 1996). The 
main advantage of CMB anisotropies as opposed to other, more local probes, is that they 
are sensitive to the universe in the linear regime, which statistical properties can easily be 
calculated starting from ab-initio theoretical models and compared to the observations. A 
dozen or so cosmological parameters could be extracted from the observations produced by 
the future satellite and ground based experiments. There are two potential problems in this 
program. First is the somewhat uncertain amount of galactic and extragalactic foregrounds, 
which could severely limit our ability to extract cosmological signal from the data. Second 
is the degeneracies among some of the cosmological parameters, which allows only certain 
combinations to be determined accurately, but does not allow to break the degeneracies 
between them (Bond et al. 1995; Jungman et al. 1996; Zaldarriaga et al. 1996). 

It is therefore important to investigate other independent confirmations of results 
produced from CMB anisotropics and it has long been recognized that polarization in 
the microwave sky might provide such an independent test (Rees 1968, Polnarcv 1985, 
Bond & Efstathiou 1987, Crittenden, Davis & Stcinhardt 1993, Frewin, Polnarcv & Coles 
1994, Coulson, Crittenden & Turok 1994, Crittenden, Coulson & Turok 1995). Like 
temperature anisotropy, polarization probes the universe in the linear regime and so can 
provide information useful to determine cosmological parameters. It is specially important 
for determining parameters which are only weakly constrained by the CMB anisotropies 
alone, such as the epoch of reionization or the presence of tensor perturbations. In 
both of these cosmic variance is the limiting factor in our capability of extracting the 
parameters, so measuring polarization would increase the amount of information and allow 
for a more accurate determination of the parameters. Both polarization-polarization and 
polarization-temperature correlation give an independent set of power spectra, which have 
different sensitivity to different parameters, so combining them results in a much larger 
information about the underlying cosmological model and can significantly increase our 
capability of extracting useful information from the CMB measurements. 

The main disadvantage of polarization is that it is predicted to be of a rather low 
amplitude, of the order of 10% of temperature anisotropy, so measuring it represents 
an experimental challenge that has yet to be overcome. Currently there are no positive 
detections of polarization, with the best upper limits of the order of 25ijlK (WoUack et 
al. 1993). This situation needs not continue in the future, however, as the experimental 
sensitivity increases and new larger and better experiments are being planned. Moreover, 
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as shown in this paper, cosmologically induced polarization has a unique signature in 
the data that cannot be mimicked by other foregrounds and provides a clear test of the 
searched signal. An experiment at Brown (Timbie 1996) plans to reach sensitivities of a few 
IxK, which should be sufficient to detect the polarization signal if our current expectations 
are correct. Interferometer observations could be able to extend this to much larger areas 
of the sky and produce maps of polarization both in real space and in frequency space. 
Unfortunately, at present none of the planned interferometer experiments is including 
polarization, although as argued in this paper frequency space has many advantages in 
search for the unique signature of cosmological polarization. Finally, all sky satellite 
maps of polarization are also planned by MAP and COBRAS/SAMBA sateUites and may 
eventually provide us with very high accuracy maps of polarization pattern in the sky. 

The outhne of this paper is the following. In §2 statistical properties of polarization 
parameters are derived in the small scale limit. Special nature of polarization induced by 
scalar perturbations allows one to use statistical methods developed in the context of weak 
lensing (Kaiser 1992). In §3 various 2-point estimators are presented both in Fourier and 
real space, together with their variances and covariances. §4 is devoted to the foregrounds 
and possible methods of their elimination. Polarization induced by tensor modes is 
discussed in §5 and the difference between the two types of perturbations is highlighted. In 
§6 theoretical predictions for polarization are computed for a variety of cosmological models 
and the ability to extract them with the satellite missions is discussed. This is followed by 
conclusions in §7. 



2. Polarization in the small scale limit 

In this section we derive the small-scale limit of temperature and polarization 
anisotropics. This limit is of special interest, because one can replace the general spherical 
expansion with the Fourier expansion and the expressions simplify considerably. The 
analysis in this section will be restricted to polarization generated by scalar perturbations. 
Tensor perturbations are discussed in §5. 

In general one can describe CMB polarization as a 2x2 temperature perturbation 
tensor T^j. Stokes parameters Q and U (we will ignore V in the following, since it cannot be 
generated through Thomson scattering) are defined as Q = T^x — Tyy and U — 2Txy — 2Tyx, 
while temperature anisotropy is just its trace, T — T^x + Tyy. The components are defined 
with respect to a fixed coordinate system (x, y) perpendicular to the photon direction n. 
The Stokes parameter Q is positive if temperature perturbation is larger along the x relative 
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to the y axis, while U parameter is positive if perturbation is larger along the upper right 
diagonal relative to the upper left diagonal. 

Equations of radiative transfer for polarization and temperature anisotropy simplify 
in Fourier space if we work in a frame with the axis defined parallel and perpendicular to 
d/d(3, where (3 is the angle between the wavevector k and photon direction n. Azimuthal 
symmetry guarantees that the Thomson scattering preserves the diagonal form of Jjj, just 
like in the case of plane-parallel atmospheres (Chandrasekhar 1960, Kaiser 1983). In terms 
of the Stokes parameters only Q is excited in this frame. This puts a restriction on the 
general form of polarization and as we show below it can be used to separate cosmologically 
induced polarization from other sources and systematic effects. Although only Q is present 
in this frame, when we rotate polarization by an azimuthal angle —0^^ to a fixed frame in 
the sky we generate both Q and U . At the observers position the expressions for T, Q and 
U are given by (Bond & Efstathiou 1987, Kosowsky 1996), 

T{n) = J d^kArik^n) 

Q{n) = J d'^kAp{k,n)cos{2(f)^^^) 

U{n) = y'c/3^Ap(fc,n)sin(20^^.), (1) 

where AT,p{k,n) are the Fourier components of temperature and polarization distribution 
function integrated over the momentum and n is the direction of observation on the sky. 
The expression for rotation angle depends both on k and n, but if we restrict our 
attention to the directions n around the pole then it can be approximated with —(pk, where 
(pk is the azimuthal angle of vector k. This approximation breaks down for wavemodes k 
close to the pole {z direction), but for sufficiently small scales the contribution from these 
modes to the total power becomes negligible. 

The solution for Aj-^pik, n) can be written as an integral over the sources along the 
line of sight (Seljak & Zaldarriaga 1996) 

At,p(^, n) = r dr/^-^^^-^'^STAk, r), (2) 
Jo 

where ST,p{k,T) are the source functions for temperature and polarization and can be 
expressed in terms of metric, baryon and photon perturbations (see Seljak & Zaldarriaga 
1996 for their explicit expressions). By combining equations |T] and ^ we can write the 
complete solutions for T{n), Q{n) and U{n). The only term that depends on the direction n 
in equation ^ is the exponential. The expressions for polarization can therefore be simplified 
by noting that for the directions n near the pole cos(20^) and sin(20^) can be written in 
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terms of second derivatives with respect to 6', the 2-dimensional projection of n onto fixed 
(x, y) plane perpendicular to the pole. This leads to 

Q{n) = DQ{n) J d'kAp{k,n), DQ{n) = {86,06, - d6yd6y)Vf 

U{n) = Duin) j SkApik.n), Du{n) = 2d6,d6yVf, (3) 

where Vg^ is the inverse 2-d Laplacian with respect to 6 and 6,, 6y are the components of 
6 in the fixed basis on the sky. Because Aj-^pi{k, n) depend only on the angle between the 
two vectors one can expand them in Legendre series, 

AT,p(fc,n) = E(2^ + l){-^Y^T,pmPl{^A. (4) 
I 

where fi = k ■ n. The rms values for AT^pi{k) can be obtained by solving the Boltzmann 
equation in differential form or the integral solution itself (Bond & Efstathiou 1987, Seljak 
& Zaldarriaga 1996). 

Each of the observable quantities T(ra), Q{n) and U{n) can be expanded on a sphere 
into spherical harmonics or their derivatives, 

= J2^Timyim{n) 

Im 

Qi^) = ^apimDQ{n)Yim{n) 

Im 

U{n) = ^apimDu{n)Yim{n) 

Im 

aT,pim = M-^y J d'kY;^{k)ATMk)- (5) 
The statistical properties of aT,pim coefficients follow from equation | above, 

i'^Tlm'^Tl'm') = ^W^mm' j d^kA3pi{k) = 5ll/5mm'CTl 



{a*pim'^PVm') — Slir6mm' J d^k/\!^pi{k) = Sll>6mm'Cpi (6) 

The cross-correlation between temperature and polarization is given by 

.ap,^') = 5.'5w / d'kAr^Ap.ik) ^ 6a^6^m'Cc, (7) 



Because we are only interested in n near the pole one can approximate the sphere 
locally as a plane, in which case instead of spherical decomposition we may use plane wave 
expansion. In this limit we replace Y^im'^PimXimij^) with / d'^lP{l)e^^'^ (and analogously for 
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temperature anisotropyQ) . Differential operators Dqln) and Du{n) acting on e*''^ become 
simple again and bring out cos(20j') and sin(20;-), respectively, where 0^- is the direction 
angle of 2-dimensional vector / with amplitude I. This leads to 

r(^) = {2ti)-^ j (fl/-^Til) 

Q0) = {271)-^ j Sl/-^P{l)cos{2(t)f) 

U0) = {271)-^ j dH/'P{l)sm{2(P^). (8) 

T(/) and P{1) are the Fourier components of temperature anisotropy and polarization in / 
space and have the statistical properties, 

{T{l)T\l')) = {2ixfCTi6o{l-l') 
{P{1)P*{P)) = {27t)^Cpi6d{1-1') 

{Tir)P*{n) = {27: fccMi-n, (9) 

where — /') is the Dirac S function as opposed to the Kronecker 6 in the discrete case 
and Ci is assumed to be a continuous function obtained by interpolation from the discrete 
spectrum defined in equation ^ 

To generate a map of temperature anisotropy and polarization one proceeds in the 
following way. For each pair of vectors /, — / on a discrete mesh one diagonalizes the 
correlation matrix Mn = Cti, M22 = Cpi, Mu = M21 = Ccu where I is the amplitude of 
vector /. One then generates from a normalized gaussian distribution two pairs of random 
numbers and multiplies them with the amplitudes given by the square root of the correlation 
matrix eigenvalues. Rotating this vector pair back to the original frame gives a realization 
of T{1) and P{1) (and their complex conjugates corresponding to — /), from which follow 
Q{1) = P{1) cos(20;-) and U{1) = P{1) sin(20^-). Fourier transform of T, Q and U back into 
the real space gives a map of these quantities in the sky in the small-scale limit, with the 
correct auto and cross correlations among all the quantities. Note that this differs from the 
prescription given in Bond & Efstathiou 1987. 



3. 2-point estimators 



"'^A somewhat better correspondence between small scale and large scale expressions is achieved if one 
uses i + 1/2 as the amplitude of a wavevector that corresponds to Ci (Bond 1996). 
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The 2-point correlations can be calculated either in angular or in frequency (Fourier) 
space. While the two are related via a Fourier transform, there are certain advantages to 
the analysis performed in frequency space. In the first subsection we explore this approach, 
while the correlation function approach and the comparison between the two are explored 
in the next subsection. 



3.1. Power spectrum analysis 



From a map of T(^), Q{9) and U{ff) we can obtain their analogs in frequency space 
by Fourier transform, X(l) — J dFOe'^^'^ X{9), where X stands for T, Q or U. Using the 
expressions given in the previous section one obtains the two-point functions of these 
quantities. 



{T{l)T*{V] 
{T{l)Q*{r 
{T{l)U*{V] 

{u{i)u\r 



= {2TlfCTl5D{l-l') 

= {2-KfCpi cos^ {2(1)^)5 d{1 - I') 
= {2nfCpism\2(P^)6D{l-l') 



(10) 



We see that the correlations in polarization give rise to a very characteristic anisotropy 

— * 

pattern in I space. This arises from the fact that polarization was not generated from a 
general mechanism, but rather through a process of Thomson scattering, which cannot 

— * 

generate U and V components in the k dependent frame defined in previous section. This 
characteristic anisotropy can therefore be used to separate true signal from instrumental 
artifacts and foregrounds, which is discussed in more detail in next section. Alternatively, if 
foregrounds can be kept under control one can use the characteristic anisotropy to separate 
scalar induced polarization from the one induced by vector or tensor modes, which do not 
obey the same anisotropy pattern (§5). 

To estimate the sensitivity that is possible to achieve in a measurement of polarization 
power spectrum let us assume the measurements arc given on a square grid of pixels with 
a total solid angle fl. The antenna beam smearing will be described with b{l). In the case 
of single dish observations with a gaussian beam this is given by b{l) = e'^'^t/^^ where is 
the gaussian size of the beam. In the case of interferometers b{l) is either 1 or 0, depending 
on whether the particular frequency is observed or not. We will assume that different 
measurements are uncorrelated by taking mesh spacing large enough to ignore correlations 
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induced by finite window (Hobson & Magueijo 1996). In the case of single dish experiments 
each pixel in real space has a noise contribiition with rms noise amplitudes a^. up for 
temperature and both components of polarization, respectively (we assume for simplicity 
that Q and U are being measured equal amounts of time). We will also assume that noise 
is uncorrelated between different pixels and between different polarization components Q 
and U. This is only the simplest possible choice and more complicated noise correlations 
arise if all the components are obtained from a single set of observations. In the case of 
Brown polarization experiment (Timbie 1996) the polarization measurements will be made 
by rotating the antenna axis by 45°, each time measuring directly the difference between 
the two orthogonal polarizations. This would therefore provide a direct measurement of Q 
and U components with no noise mixing between them. In the case of interferometers each 

— * 

pixel in / space is measured directly and the noise is uncorrelated between individual pixels 
in frequency space. Following Knox 1995 we will introduce pixel independent measure of 
noise p = flapp/N for single dish experiments and for interferometers 

(see Hobson & Magueijo 1996 for expressions that relate crr,p to the receiver sensitivity). 
For cross-correlation between temperature and polarization there are two simple cases 
to be considered. In one the cross-correlation is being made with two different maps, in 
which case noise is uncorrelated. In the second case both temperature and polarization 
are obtained from the same experiment by adding and differenciating the two polarization 
states. In this case noise in temperature and polarization are related via — (t%/2. In 
both cases noise in temperature is uncorrelated with noise in polarization components, so 
the final expressions are identical. 

The first step is to construct a discrete Fourier transform of the map, 
X{1) = N-^Y.e~^'^'^X{e). Using the observed quantities f(l), Q{1) and Uil) we can 
form power spectrum estimates. The simplest one is for temperature anisotropy, which for 
single dish observations is given by 



where the term h~'^{l) accounts for beam smearing and A^^ is the number of independent 
modes around I. The expression for interferometers is the same without the beam smearing 
term and summing only over the observed modes in frequency space. We will only present 
expressions for single dish here, as the modification for interferometers is obvious. Each 



If we average over Ni modes the variance is reduced by A^^ . There are /^f2(iln//27r modes 






r'. The 

(/) + tVj 
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of amplitude / in an interval d\nl, which leads to the variance 



cov(c; 



Tl) 



^2 47r 



(12) 



where Gov (XX') = ((X — (X))(X' — (X'))). If there is more than one field the variance 
decreases inversely proportional to the square root of number of fields (Hobson & Magueijo 
1996). In the limit of large sky coverage it reduces to the expressions given by Knox 1995 
and Jungman et al. 1996. 

In the case of polarization there are several estimators that one can form. The simplest 
one is given by 



C 



PI 



I 



-1 
p 



(13) 



The single mode pair variance for this estimator is {Clih\l) + 2Cpib\l)wp^ + 4wp^y/\ If 
noise is dominant contributor to the variance, as it is likely to be the case for polarization 
given the small overall amplitude of the signal, then this estimator is far from optimal. The 
optimal estimator is 



C 



PI 



n 
'n, 



J2 \Q{1) cos{2(f)^) + U{1) sin(20,-3|2 - w-p^ 



I I 



b-\l), 



(14) 



which has single mode pair variance CilP'il) + Wp^ . In the limit of noise dominated variance 



this is 2 times smaller than the variance of the estimator in equation |T3|. Averaging over all 
modes in dXnl gives 



cov(c; 





"/2Qdln/ 



'\ + \wpCpih\l)\-^\ 



(15) 



Finally, for cross-correlation the optimal estimator is, 

C'c/ = ^ E + Q\l)m) cos(20,-) + {ij{S)f\l) + lJ\l)f{S)) sin(20^-)] h-\l), 

(16) 

which has a single mode pair variance \{C}jp'^{l) + (Crib'^il) + Wrp^){Cpib'^{l) + typ^))/2]^/2. 
As before, averaging over the modes reduces this variance inversely proportional to the 
square root of the number of modes and gives. 



Cov(C, 



ci) 



C, 



271 



1 + 



iCTib\l)+WT')iCpib\l) + Wp') 
Chb\l) 



;i7) 
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For a study of cosmological parameters one also needs to include the covariance 
elements between various power spectrum estimators. These are given by 

Cov(C.,C„) = clM^^, 

Att 

Coy{CciCti) = CciCTi-^^^^{l + [wTCTih\l)]-^) 

4.77- 

CoY{CciCpi) = CciCpi^^^^^{l + [wpCpih\l)]-'). (18) 



3.2. Correlation function analysis 



In this subsection we explore the correlation function analysis of CMB anisotropics. 
Taking Fourier transform of the power spectra leads to the following correlation functions, 



(T(O)T(^) 
(Q(O)Q(^") 

{Qm0) 

imQie) 



^-^CTib\i)Mie) 

ZTT 

J ^-^Cpib\l)[Jo{W) +cos{A(P) J,{W)] 
'■ Idl 

—Cpib\l)[Jo{W) -cosm)J,{W)] 

4:71 

^-^Cpih\l)sin{A(t>)UW) 
J^-^Ccib\l)cos{2<P)J2{W) 

^-^Ccib'{l)sm{2<P)MW), 
47r 



(19) 



— * 

where is the direction angle of 9 and J„(x) are the Bessel functions of order n. 
Although the characteristic anisotropy is present also in the correlation functions, it is 
more complicated, because there are actually two independent correlation functions for 
polarization. 



cue) = J^-^Cpib'{i)MW) 
idi 



cm 



An 



Cpib\i)j^{ie), 



(20) 



both of which of course depend on the same underlying power spectrum. Moreover, the 
prior Cp; > becomes a set of integral constraints in real space, which is more difficult to 
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impose on the estimators obtained from the data. All this argues for Fourier space analysis 
as the method of choice in the case of polarization. 

Expressions above agree with those derived by Coulson et al. 1994 and Kosowsky 1996, 
which for QQ and UU are also only valid in the small scale limit, but where this limit has 
not been consistently applied to all the steps, so that their final expressions look much more 
complicated than they actually are. They contain a term involving double summation over 
Ap;, Ap;/ in the (f) dependent term, which as shown in Appendix reduces to the expressions 
above if one consistently applies small scale limit to their expression. All the information 
about polarization is therefore contained in Cpi and Cci- Note that taking the Fourier 
transform of QQ or UU correlation function does not result in the power spectrum of Cp;, 
but has an additional dependent term that involves a double integral of Jq{W) and J^^il'O) 
over /' and 6. Although generally smaller than Cp/, this integral does not vanish in general, 
hence the appearance of such terms in expressions by Coulson et al. (1994) and Kosowsky 
(1996). This is of course not the optimal way to obtain the power spectrum from the 
correlation function. To obtain an estimate of the underlying power spectrum it is better 
to work in the Fourier domain directly, following the methods given in previous subsection. 



3.3. Predicting polarization from temperature maps 

One can use the measured temperature maps to predict the polarization pattern. The 
estimator is 

Q{1) = aT{l) cos(20;-.), U{1) = aT{l) sin(20;-), (21) 
where by minimizing the variance 

mi) - Q{l)f) = cos\2<Pj)[a'CTi - 2aCci + C„] (22) 

and similarly for U one finds a = Cci/Cti, so that the fractional variance in the estimator 
is 

<MzW)]!>=i_Corr(T.P)f. (23) 

The correlation coefficient is defined as Corr(T,P)i = Cci/(CtiCpi)^/^. Figure |l|d shows 
that the correlation coefficient typically ranges between —0.5 and 0.5 and so the fractional 
variance will be at best around 0.8 or so in / space and even larger than that in real space, 
where one averages over positive and negative cross-correlations in power spectrum. This is 
not very impressive in terms of predicting where to look for large polarization amplitude. 
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although in a statistical sense one is still much more likely to find a high signal at a high 
(7 peak in the temperature map than at a random point (Coulson et al. 1994). On large 
angular scales (/ < 10) the correlation coefiicient can be much larger and approaches 
unity in some models, so that the expressions in this limit would actually give a good 
correspondence between observed and predicted polarization. Unfortunately, the amplitude 
of polarization is extremely small on these scales and there is little hope to measure it in 
the near future even with the help of this "matching filter" technique. 



4. Removal of foregrounds and other systematics 

Although several galactic and extragalactic foregrounds are significant in the case of 
temperature measurements, only few of these are polarized and need to be considered for 
polarization measurements. Radiation from earth atmosphere and bremsstrahlung emission 
are not polarized at millimeter wavelengths and need not be discussed further (although 
atmospheric emission does have an effect by increasing the effective temperature of receiver 
and adding a fluctuating offset). 

Extragalactic radio sources have synchrotron radiation as the dominant emission 
mechanism and can be 20% polarized. Their contribution to the polarization signal will be 
similar as their contribution to the temperature anisotropy signal and so analysis of point 
source effects on CMB can be directly apphed to the case of polarization as well. This is 
discussed in more detail by Prancheschini et al. (1989, 1991) and Tegmark & Efstathiou 
(1996). Point source contamination depends on the observed frequency, angular scale and 
flux cut above which point sources can be identified and eliminated. Poisson distribution 
produces a white-noise spectrum and at large angular scales radio point sources do not pose 
a significant problem. For example, on angular scales above 1° their contribution is less 
than IfiK at 30 GHz (Timbie 1996), which is below the expected amplitude of the signal 
and is even less than that at higher frequencies. On smaller angular scales point sources 
become more important and more ambitious fiux cuts are needed, which limits the area 
of the sky that can be observed. While in the case of temperature anisotropy this can be 
the main limitation of an experiment (and indeed of the whole CMB field), in the case of 
polarization we do have an additional constraint that allows us to separate cosmologically 
induced signal from the foregrounds. This is discussed in more detail below. 

On large angular scales the main foreground contribution comes from our galaxy, where 
both dust and synchrotron emission can be polarized. Dust emission in the far-infrared 
is polarized up to 10% (Hildebrand et al. 1995). The contribution of dust to the lower 
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frequency channels (where most HEMT based polarization measurements are being 
planned) is small. At frequencies below 100 Ghz the most important source of polarization 
is likely to be galactic synchrotron emission. Its linear polarization can reach 70%, so 
that polarization amplitudes of 50/iK are expected around 30GHz (Cortiglioni & Spoelstra 
1995). This drops significantly at higher frequencies and only a few fiK synchrotron 
contamination in polarization is expected around lOOGHz (Timbie 1996). Nevertheless, 
this is of the same order as the expected signal, so that some further foreground rejection 
is needed. One possibility is to use only clean parts of the sky where synchrotron emission 
is low, such as at high galactic latitudes. In addition, multifrequency CMB observations 
can be used to remove the foregrounds, which in the case of only one important foreground 
with approximately known frequency dependence can be very effective (Brandt et al. 1994). 
Accuracy of l/iK can be achieved with only two frequency channels if the noise level is 
around 1/iK per pixel (Timbie 1996). Another possibility is multifrequency removal using 
the more sensitive temperature maps. This would be a useful strategy for example in 
the case of COBRAS/SAMBA satellite, where only lower frequency channels will have 
polarization capabilities, but all frequency channels could be used to determine the local 
contribution of various foregrounds. This way one could effectively remove multicomponent 
foregrounds from polarization even if only a few channels actually measured it. 

While each of the foregrounds above can in principle be removed from the data, in 
practice this may not always be possible at the levels of 1/iK and it would be useful to 
have an additional test of the presence of cosmological signal. Characteristic anisotropy 
of polarization in Fourier space provides such a test and gives a unique signature of 
polarization induced by scalar perturbations. To test whether the signal is cosmological one 
needs to compare the quantities 

E{f) = Q{1) cos{2(t)j) + UiJ) sin(2(/);-) (24) 

and 

B{T) = -Q{1) sin(20^-') + U{1) cos{2(j)^). (25) 

The first quantity contains all the polarization signal and its estimator is given in equation 
T3, while the second quantity vanishes even in the presence of cosmological polarization 



induced by scalar perturbations. This is true not just statistically, but for each Fourier 
mode individually. Most of the foregrounds should contribute on average the same amount 
to both variables. This is certainly valid for uncorrelated point sources, but even in 
the case of synchrotron and dust emission the alignment is preferentially determined by 
magnetic fields, which are not scalar in nature and will not exhibit the characteristic 
anisotropy in Fourier space. Hence the difference between E and B can be taken as a 
measure of the cosmological signal as compared to the foregrounds and/or instrumental 
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offsets. If foregrounds and not noise are expected to be the main limitation then one may 
use may subtract the power spectrum of E from the one for B in equation |1^ and no 
further foreground removal is needed. One could therefore use this technique to measure 
polarization even at frequencies below 50GHz, where foreground contribution is large, but 
could be averaged over if sufficient number of channels are being measured. Note that this 
test does not depend on the temperature anisotropics at all and can be apphed directly to 
the measurements of Q and U. To obtain a statistically significant measure of polarization 
one only needs to show that in an rms sense E is larger than B. If one is noise and not 
foreground limited then it is more advantageous to use the optimal estimator in equation 



14, because subtracting the power spectrum of B from the optimal estimator for E leads 



to an increase in noise by 2^/^. In practice the actual analysis will depend on the level 
of foregrounds and other systematics (e.g. sidelobe pickup) relative to noise and different 
estimators must be tested for consistency, but it is important to note that in the case of 
polarization we have a possibility to use a combination of Stokes parameters in which 
foregrounds can be separated from the cosmological signal, something that cannot be 
achieved in the temperature measurements alone. 



5. Tensor polarization 



Discussion in §4 applies only if polarization is produced by scalar perturbations. While 
this is certainly a valid approximation on small angular scales (/ > 100), on larger scales one 
may be able to detect polarization from non-scalar perturbations, induced either by vector 
or tensor modes. The latter are of particular interest, because they are expected to be 
present in several infiationary based models, although only on large angular scales and with 
rather small amplitudes (Crittenden et al. 1993). Defect models also predict production 
of both tensor and vector perturbations. In the case of tensor perturbations temperature 
anisotropy and the two Stokes parameters can be written in the small scale limit as (Bond 
1996, Kosowsky 1996), 

T(^)(^) = /rf'^(l-/i')Af+)(fc,/i) 

Q^^^ (n) = I (fk [(l + /i^) Ag'^^ (fc, /i) cos(2(/)^^) - 2/iAg'^) (fc, /i) sin(20^^)' 

(n) = j d'k [(l + /i') Ag^+^ {k, fi) sin(20^) + 2/iAg^^^ {k, /i) cos(20^)] , (26) 

where Ap'~^\k, jj,) and Ap'^\k, jj,) are the two independent polarization states of a gravity 
wave with equal rms amplitudes. For convenience we defined the orientation of the two 
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polarization states in the plane perpendicular to k so that local x direction points in the 
direction of the pole, in which case only A^- {k,fi) contributes to the temperature in the 
small scale limit. Expectation values for A^p~^\k,jj,) and A^p^\k, /j,) can be calculated just 
like in the scalar case by expanding them into a Legendre series and solving a system of 
Boltzmann equations (Crittenden et al. 1993) or the integral solution (Zaldarriaga & Seljak 
1996). Because of additional n terms in equation ^ a more convenient set of variables is 
obtained by eliminating the explicit fi dependence (Kosowsky 1996), 



21 + 1 
1 

2/ + 1 



(/ + l)Ag,ti + ^A(^^^ 



pi-i 



(/ + l)(/ + 2) (T,) ^ ^ 6/^ + 9/2_/_2 ^(^,) , {1-1)1 ^ 

iApj,r,-\-^-7-: ———-—lApj i — —LA 



21 + 3 -^''+2 ^ - - 1)(2/ + 3) "^'^ 2/ - 1 



(27) 



where e stands for + and x and all the variables explicitly depend on k. 



One can now follow the same steps as in the case of scalar perturbations, which 
transform the angle 0^ into 0^-. The variables E^'^'^ and B^'^^ defined in equation ^ (where 
superscript T indicates that these are produced by tensor modes) are independent of the 
azimuthal angle and the two tensor components decouple, so that Ap {k,fj,) contributes 
only to and A^p^\k, jj) contributes only to B^'^\ Their power spectra can be expressed 
in terms of integrals over B}, Bf as, 

{B^'^\l)B^'^>{V)) = {2Txf5{l-l')jd'k\B]Wk) 
(E(^)(/>(^)*(/0) = i2n f 6(1- f) J d'k\Bf\\k) 
(E(^)(r)5(^)*(/0) = 0. (28) 

The cross correlation term vanishes because the two tensor polarization states are 
independent. The variable -B^^-* does not vanish in the case of tensor perturbations and its 
power spectrum differs from the power spectrum of E^'^\ Note that different combinations 
of Q and U will result in different power spectra, which can always be expressed in terms of 
the two defined above. Just like the cross term between E^'^'> and B^'^'> vanishes so does the 
cross correlation term between T*^^) and B^^\ There is only one power spectrum present in 
the case of temperature-polarization cross correlation, 

{T^^\l)B^^^*{U)) = 

(T(^)(0^^^^*(^0) = {27ry6{l-l') f d^kA^^\k)Bf{k). (29) 



Detailed calculations of these spectra have been presented elsewhere (Seljak & Zaldarriaga 
1996b). 
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6. Model predictions 



Instead of the power spectrum Ci we will use the quantity /(/ + l)C;/27r, which gives 
the contribution to the variance per logarithmic interval of /. This is a familiar quantity in 
the case of temperature anisotropies, where its broad band average is approximately fiat 
up to the damping scale. Predictions for /(/ + l)Cpi/2'K and /(/ + l)Cci/27r are given in 
figures |l]a,b for a variety of cosmological models. For comparison we also plot the usual 
/(/ + l)CT'i/27r in figure |l|c, as well as the correlation coefficient in figure |I|d. All the 
model predictions have been computed using CMBFAST package by Seljak & Zaldarriaga 
(1996). One can see that typically the models predict very little polarization on large 
angular scales, below / ~ 200. On smaller angular scales most of the models predict 
polarization at the level of 5-10/iK. There are several characteristic features of interest in 
this regime. The most important one is that the acoustic peaks are narrower than the 
corresponding ones in temperature anisotropy. One can understand this with the help of 
tight coupling approximation (Hu & Sugiyama 1995; Seljak 1994; Zaldarriaga & Harari 
1995). The dominant source of temperature anisotropy are intrinsic photon anisotropy 
(Ao) and velocity (Ai). Both terms oscillate, but are out of phase with each other. This 
means that they partially cancel each other and oscillations in the temperature anisotropy 
are less pronounced than they would be if only one term were contributing. The dominant 
source of polarization is photon second moment A2, so the oscillations are more pronounced 
than in the case of temperature anisotropy. These oscillations are even more pronounced 
in the case of temperature-polarization cross-correlation, which can be either positive or 
negative. Another characteristic of polarization is that it is not sensitive to the integrated 
Sachs- Wolfe term. This term is responsible for increase in the temperature anisotropies 
at low /, as in models with cosmological constant, curvature or if recombination occurs 
close to the matter-radiation equality, during which gravitational potential is changing with 
time. Another class of such models are topological defect models, where small I spectrum is 
dominated by late time integrated Sachs- Wolfe effect. A measurement of polarization at a 
few microkelvin level will put a significant constraint on such models. 

One of the parameters that are of special importance for polarization is the optical 
depth to Thomson scattering. As photons propagate through the universe they scatter off 
free electrons, which were ionized by UV light either from an early generation of stars or 
from quasars. Current limits give that the universe was mostly ionized up to the redshift 
of 5, which results in optical depth of the order of 1% in standard CDM universe and 
somewhat larger in open or high baryon models. It is likely that the reionization did not 
occur much earlier so that the optical depth would exceed unity, because then it would 
suppress CMB anisotropies on small angular scales, in contrast with recent observational 
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data (Netterfield et al. 1995, Scott et al. 1996). The exact epoch of reionization remains 
however unknown and its determination would provide an important constraint to the 
models of galaxy formation. Temperature anisotropics alone will not be able to constrain 
this epoch significantly, because even in the most optimistic scenario sensitivity to optical 
depth T is around 10-20% (Jungman et al. 1996). This is because reionization is degenerate 
with the amplitude of fluctuations, which can be only broken at low I, where cosmic variance 
is large. Polarization can help here both because reionization introduces new structure 
and also simply because the cosmic variance is reduced as more independent realizations 
are observed. As shown in figure |I|, the effect of reionization is to increase somewhat the 
amplitude of polarization at low /, but not by much and the amplitude still remains below 
1-2^K. The effect is better seen in the cross-correlation spectrum and in the corresponding 
correlation coefficient. The latter clearly displays the rich structure at low 1 that allows 
one to determine epoch of reionization and the integrated optical depth (Zaldarriaga 1996). 
On smaller angular scales polarization amplitude decreases with optical depth just like the 
temperature anisotropy and the ratio of the two remains constant (Bond & Efstathiou 
1987). 

Figure presents a more quantitative estimate of sensitivity in the case of satellites, 
assuming standard CDM as the underlying model. The middle curve shown is the 
underlying theoretical model while the two curves above and below show the one standard 
deviation of the reconstructed spectrum from the true model. The variances were obtained 



using equations [T^ and |T^, assuming 50% sky coverage {^l = 2tt) and d\nl = 0.2. We 



— 1/2 

adopted 3 different noise characteristics. The most optimistic possibility is Wp =0.1 /xK 
in 0.2° beam, which could easily be achieved by COBRAS/SAMBA satellite with their 

— 1/2 

bolometric detectors. The intermediate model assumes Wp = 0.2 /iK and 0.3° beam, 
which could be feasible with the MAP satellite by combining their most sensitive channels. 

— 1/2 

The third model is the most conservative one and assumes Wp = 0.3 /xK and 0.3°. All 
of the sensitivities assume 1 year of observation and longer observation periods would 
reduce the noise accordingly. From figure |^ we see that only the most optimistic model is 
capable of constraining the polarization power spectrum significantly. On the other hand, 
for cross-correlation spectrum the situation is much better and all of the assumed models 
will give some positive detection, although of course with lower noise levels one will be 
able to extend this to much smaller angular scales. This difference between polarization 
and cross-correlation is to be expected, because noise in temperature is lower and because 
cross correlation power spectrum has more power than the polarization power spectrum 
itself. Although more detailed calculations are needed to estimate the sensitivity for actual 
satellites (Zaldarriaga et al. 1996), it is clear that reducing the noise by a factor of 2 leads 
to a significant improvement in the sensitivity to polarization. 



7. Conclusions 



Polarization in cosmic microwave background has the promise to become the new 
testing ground of processes in the early universe quite independent of the measurements 
of temperature anisotropics. The spectrum of polarization, while lower in signal than the 
temperature anisotropics, can be more sensitive to certain parameters such as reionization 
or gravity waves. Even for the determination of the standard parameters polarization 
provides some advantages, for example, acoustic oscillation peaks are much more prominent 
and so can be more easily detected. In addition to the spectrum of polarization one can 
also determine the spectrum of temperature-polarization cross-correlation. Two additional 
spectra can help to break some of the degeneracies present in the estimation of cosmological 
parameters. This is specially important for those parameters whose precision is limited by 
cosmic variance, because polarization provides additional independent realizations of initial 
conditions in the universe. 

Even though measurements in polarization are not likely to achieve the same level 
of precision as in the case of temperature anisotropics, polarization still has some 
advantages that may prove crucial if the amplitude and complexity of galactic foregrounds 
and extragalactic point sources have been significantly underestimated. Multifrequency 
subtraction is simpler in the case of polarization, becaTise fewer foregrounds are polarized 
and need to be modeled. More importantly, polarization induced by scalar perturbations 
has a unique signature in Fourier space and by exploiting this one may separate cosmological 
signal from other sources of polarization. Although the analysis in this paper has been 
limited to small scales it has recently been shown that the same property is also valid in 
the more general all-sky analysis (Zaldarriaga & Seljak 1996; Kamionkowski, Kosowsky & 
Stebbins 1996). This signature would be especially important if the level of foregrounds is 
significantly larger than expected. While one would not be able to remove the foregrounds 
from the temperature maps, a combination of Q and U Stokes parameters would allow one 
to subtract statistically the effects of foregrounds the case of polarization. This technique 
would be feasible both for interferometer measurements or for measurements with a large 
coverage of the sky such as the forthcoming satellite missions. The sensitivity of satellites 
will be sufficient for an unambiguous determination of polarization and indeed polarization 
may prove to be crucial to break some of the degeneracies present in the parameter 
reconstruction from the temperature anisotropics alone (Zaldarriaga et al. 1996). Finally, if 
foregrounds can be controlled, then a unique signature of tensor (or vector) perturbations 
could be directly observed, although this would require an exquisite understanding of noise 
properties, systematics and foregrounds at the level of 0.5/.j,K . Given the unique nature 
of information present in the microwave background and its simple linear depence on the 
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underlying cosmological parameters it is important to explore it at its maximum, which 
certainly includes polarization as one of its main components. 

I would like to thank Douglas Scott, Arthur Kosowsky, David Spergel, Martin White 
and specially Matias Zaldarriaga for useful discussions. 



A. Appendix 



In this Appendix we show how the small scale limit of correlation functions QQ and 
UU (equations |T9|) follows from the expression derived by Coulson et al. 1994. We start 
from their expression (correcting for missing factors of 1/2), 



2/ -I- 1 

{Qmm = 

(^-4)! 
(^ + 4)!i7 



CpiPi{cos9) + ^^^^CppiP^ (cost 



C 



ppi 



E(2/' + 1) / edkAp,Api4i, (Al) 



where Pi is the Legendre polynomial and P^ the associated Legendre function of 4th order. 
Coefficients afi, = ^\dxPi{x)P^{x) have a closed form expression (Coulson 1994) and in 
the large / limit they peak at / = More importantly, is a rapidly oscillating function of 
k and for / ^ V the integral over k leads to almost complete cancellation of / k'^dkApi^Api. 
Thus in this limit one can write 

2/(/-l)(/-2)(/-3)(/-4)! ^ 
Cppi ^^^-^yj Cpi, {A2) 

where we used the / = I' closed form for of af;, (Coulson 1994). Furthermore, in the / oo 
limit Legendre functions can be written as P™(cos0) = J„i{l9)l~^{l + m)!/(/ — m)!, where 
Jm{x) is the Bessel function of order m. Combining all the expressions leads to {Q{0)Q{0)) 
correlation function given in equation |l^. Other correlation functions in equation |19| can be 
derived from the expressions in Coulson et al. 1994 using similar manipulations. Numerical 
results presented in Zaldarriaga & Seljak 1996 show that small scale approximation derived 
here is an excellent approximation everywhere except at very small / (/:$ 20). 
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Fig. 1. — Power spectra of polarization (a), temperature-polarization cross correlation (b), 
temperature (c) and correlation coefficient (d). The models are standard CDM (solid curve), 
open CDM (dotted curve) and reionized standard CDM with optical depth of 0.2 (dashed 
curve). 
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Fig. 2. — Variance in polarization (a) and cross-correlation power spectrum (b) for a satellite 
with noise characteristic wp = wt/^. We assumed 50% sky coverage and 0.3° degree beam 
(0.2° in the most optimistic case). The spectra were averaged over a 20% band in / and the 
bands shown are one standard deviation above and below the underlying model, taken to be 
COBE normalized standard CDM. 



